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Abstract. A Generalized Kinetic Theory [1] was proposed in order to have the 
possibility to treat particles which obey a very general statistics. By adopting the same 
approach, we generalize here the Kinetic Theory of electrons and phonons. Equilibrium 
solutions and their stability are investigated. 

1. Introduction 

Very recently [1] , a Generalized Kinetic Theory (GKT) has been proposed by Rossani 
and Kaniadakis, in order to have the possibility to treat, at a kinetic level, particles which 
obey a very general statistics. 

The quasi-classical Boltzmann equation introduced in ref. [1] is a generalization of 
the Uehling-Uhlenbeck equation [2], which is intended for bosons and fermions only. In 
ref. [1] we have shown that our generalized Uehling-Uhlenbeck equation (GUUE) assures 
particle, momentum and energy conservation. Equilibrium, its uniqueness and stability 
(via an H-theorem) have been investigated. 

In ref. [3] Koponen points out that fractal or inverse power law distributions of phonon 
excitations are of interest in modeling various meaningful situations in solid state physics. 
Moreover, he feels that until recently there has been little guidance on how to generalize 
the kinetic theory of electrons and phonons obeying non Gibbsian statistics. 

Here we propose a Generalized Kinetic Theory for Electrons and Phonons (GKTEP), 
by following the same ideas which lead to the GUUE. In order to keep our GKTEP as 
general as possible, we introduce modified collision terms not only for phonons, but also 
for electrons, so that an application is allowed not only to electrons, but also to other 
particles (obeying a general statistics) which interact with a crystal lattice. 

The paper is organized as follows. In section 2 the Bloch-Boltzmann-Peierls [4] equa- 
tions are briefly recalled, and our generalization is introduced. 

In section 3 and 4 we study the equilibrium solutions to the GKTEP equations and 
their stability is investigated, via an H theorem. 

2. Generalized kinetic equations for electrons and phonons 

The most complete description, at a mesoscopic level, of a system of electrons and 
phonons, is based on the Bloch-Boltzmann-Peierls (BBP) equations [4]. 

Consider two populations: electrons (e), whose number is conserved, obey Fermi-Dirac 
statistics; phonons (p), whose number is not conserved, obey Bose-Einstein statistics. Of 
course, one should deal, in principle, with p-p, p-e, and e-e interactions. A first important 
assumption is the following: electrons (distribution function n) are considered as a rarefied 
gas in a "sea" of phonons (distribution function A^). This means that e-e interactions 
can be neglected. In e-p interactions only electrons and energy are conserved. In p-p 
interactions only energy is conserved. 
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By adopting the notation of ref. [5], let Ng — Ng{\s.,yi,t) be the distribution function 

of phonons (quasi-momentum k, energy ujg{k.)) of type g (i.e. branch g of the phonon 
spectrum) and n = ?i(p,x, t) the distribution function of electrons (quasi-momentum p, 
energy e(p)). The distribution functions Ng and n are normalized so that the thermal 
energy density Ep of the crystal is given by 



while the concentration N and the energy density Ef, of the electron gas are given, respec- 
tively by 

= J ^(P,x,t)2dp, £^e(x,t) = ^ j e(p)n(p,x,t)2dp, 

where the factor 2 inside these integrals accounts for degeneracy. 
The BBP equation for phonons reads 

(ug = dojg/dli), where 

(^) = I ^-(V2) E ^ - - ^9.)[Ng{Ng, + l){Ng, + 1) 

^ ^ VP J g,g^ 

- {Ng + mgiNg,] + J2 ^ppO^^ ki ^ ks)S{u;g, -Ug- Ug,)[{Ng + l){Ng, + l)Ng, 
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NgNg,{Ng, + 1)]} 



dNg 



J 2wpe{p p', k)[np(l - np/)(l + Ng) - np/(l - np)Ng]5{ep^ +^9- ^p) 



dp 

8^' 



where Up — n(p), ep = e(p), Wpp and Wpe are transition probabilities. 

The first term in the braces corresponds to the direct and reverse processes 

^ («/i,ki) + (^2,k2), 

where 

k2 = k - ki - b (2) 
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and b is a vector of the reciprocal lattice. 

The second term in braces correspond to the processes (direct and reverse) 

(53, ks) ^ (^,k) + (^i,ki), 

where 

kg = k + ki + b. (3) 
The BBP equation for electrons reads 

(e-e interactions have been neglected), where — e is the electron charge and 

(^) = E/ {«^ep(p',k- p)[np,(l -np)iV, 
- np(l - np')(l + Ng)]6{ep + Sp' - oug) 



+ Wepip' p,k)[np'(l - np)(l + Ng) - np(l - np')Ng]S{€p + u)g - €p')} 



where Up = n(p), v = dep/dp, and the transition probabilities tUep's obey the following 
relationships: 

Wpe{p P', k) = Wep{p p', k) = Wepip', k p). 

The first term corresponds to processes with emission of a phonon having quasi- 
momcntum k by an electron having a given quasi-momentum p, and reverse processes 
with absorption of a phonon k by electrons p' with return to the quasi-momentum p: 

p = p' + k + b. (5) 

The second term corresponds to processes with absorption of a phonon by an electron 
p and the reverse processes of its emission by electrons p': 

p + k = p' + b. (6) 

We shall generalize [1] the expressions of 

dNg\ fdNg\ fdrip 



dt Jpp V dt J^^ V 5^ / ep 



by introducing the following substitutions: 



1 + iV, ^ m{Ng), Ng ^{Ng) 

l-Up ^ i>{np), Up (p{np), 



where ^, and (f, are, respectively, non negative functions of A^^ and Up, obeying the 
conditions 

^r(O) = ^/;(0) = 1, $(0) = ^(0) = 0. 

We assume that ^{Ng)/'if{Ng) and ip{np)/i/){np) are monotonically increasing functions, 
respectively, of Ng and Up. This is trivially true for bosons and fermions. In general, this 
assumption is justified a posteriori, since it assures uniqueness and stability of equilibrium. 

The functionals ip, (p and $ replace the Pauli bloking and boson enhancement 
factors in the phase space, respectively. Alternatively, one could say that generalized 
statistics, defined by means of the couples "if, $ and i/j, (f, are introduced for phonons and 
electrons, respectively. 

We have then 



-(1/2) J {J^WppO^iM - k)5(a;, - a;,, -a;,J[$(iV,)*(iV,J*(iV,J- 



dt , 
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^{Ng)^NgJ^Ng,)] + ^«;^^(k,kl ^ k^) S g, - U g - U g,) {N 9i) m ' 
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m9)m9i)m9s)]}^: 



dNg 



I 



dt 

2wpe{p,p' k)[(p{np)'ip{np')'^{Ng) - (p{np')'ip{np)^{Ng)]5{ep' + Ug - ep) 



pe 

dp 



dup 
~dt 



= J]/ {wep{p\^^p)[v{np')il}{np)^{Ng) 



ep g 

- v{np)'4^{np')'^{Ng)]5{ep + ep/ - ujg) 

dk 

+ Wep{p' p,k)[(^(np/)^(np)*(iVg) - ip{np)i){np,)^{Ng)]5{ep + ujg- v)}^- 
3. Equilibrium 

Even though the generalized equations are certainly more complicated that the original 
ones, many of the usual techniques of kinetic theory are still applicable in the study of 
equilibria and their stability. In particular, we will prove an H theorem for the present 
problem, whose connections with the entropy law are easily established. 

In the space homogeneous and forceless case, equilibrium is defined by 

_ f^^p \ =0. (7) 



9t J coll \ 9t / CP 

In order to study equilibrium and its stability, it is useful to introduce the following func 
tional: 



^ I ^ dt J^^,, ^{Ng) J \dt VK) 



which can be put in the following form: 
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+ ^9.^92^3) In I'^l'^l'' dk,dk2 + 2J2 I I wep{p^p',k)x (9) 



^91^92^93 g 




where 



X (-<^V$g + <^V'*5)<^(ep' +^9- ep) In ^^^^ dkdp < 0, 



=$(iV,(k,)), vEr^, = *(7V,(k,)), 

if' = (p{np'), ip' = ip{np>). 



Proposition 1. 

Condition (7) is equivalent to the following couple of equations: 

vE'(iV,)$(iV,J$(iV,J = $(iV,)$(iV,Jvl>(iV,J Vk, ki (10) 
cp{np),p{np,)^{Ng)^if{np>),p{npMNg) V p, k. (11) 

Proof. First of all we observe that (10), (11) (7). On the other hand, from (8) 

we have (7) =^ V = 0. Now, since in eq. (9) both the integrands are never positive, one 
has P = =^ (10), (11) • 

Condition (10) shows that ln($*/\E'*) is a coUisional invariant for phonons, that is 

$(iV„*) 

In ] '^ i ^ -uJT (12) 

(hereinafter * means "at equilibrium"), where T is the absolute temperature of the whole 
system (electrons plus phonons). 

From (11) and (12), taking into account that ep = ep'+a;^, we find that, at equilibrium, 
ln[<^(np)/'0(np)] + e^/T is a coUisional invariant for electrons, that is 

ln^ = (/x-6p)/T, (13) 

where ji is the chemical potential of the electron gas. 

Observe that, due to the monotonicity of both (^/^ and ^p/ip, equations (12) and (13) 
give unique solutions for N* and n*, respectively. 

We would like to stress that, with a proper choice of the ratios ^»/^' and v'/V', it is 
possible to reproduce the non-standard quantum statistics proposed in the last years and 
compare them from the point of view of the GKTEP, but this will be matter of a separate 
future paper [6]. 
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4. Stability 

In order to study the stability of such equihbrium solutions, let us introduce the 
following functional: 



where 



d-HpiN,) _ $(iV,) dn^jn^) ^ .v(^p) 



dNg ^{Ng)' dup iPirip)' 

Observe that, since and (f/ip have been assumed to be monotonically increasing. Hp 
and He are convex functions of Ng and Up, respectively. 

We can now prove the following H theorem: 

Proposition 2. 

H is a Lyapounov functional for the present problem. 

Proof. First of all, it is easily verified that H = V < 0. Then we observe that 

where Hpg — Hp{Ng) and Hep = He{np), due to electron and energy conservation. 
Now we can write 



H-H* = J2j npgdk + J Hepdp, 



where 

'l-a — 'l-a 



where (/>q, = Ng for a = pg and = Up for a = ep. 

Due to the convexity of Ha, we have Ha > and finally H — H* > • 
By definining entropy as 

o _ c _i_ c _ 

^ - ^P + ^e- ^^^3 , 

these results are easily interpreted on a physical ground. In fact, H < and H — H* > 
simply mean that entropy is always non decreasing and it attains a maximum at equilib- 
rium. 

Observe that this definition of S is consistent with the definitions of T and fi we 
already gave. In fact, at equilibrium, the following thermodynamical relationships are 
recovered: 

where E = Ep + Ee- 
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